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The mechanism at the base of phase noise generation in electrical oscillators is reexamined from 
first principles. The well known Lorentzian spectral power distribution is obtained, together with a 
clearcut expression for the line-width parameter. The mechanism of suppression of the amplitude 
fluctuations and its effects are discussed and the true role of the figure of merit Qo of the resonator 
is restated. The up conversion of low frequency components from a non-white noise source is also 
considered. A number of simple numerical experiments is presented to illustrate and clarify the 
mathematical results. 


INTRODUCTION 

Signal generation is half the world of electronics, the 
electronic oscillator is ubiquitous, it has a part in all con¬ 
sumer, professional and scientific appliances and most of¬ 
ten the quality of the generated signal is crucial to the 
overall performances. 

After the problems of accuracy, stability, adjustability 
and resolution, one of the main concern in today signal 
generation is spectral purity, either harmonic content or 
non-harmonic wide and narrow band noise. 

The fact that a sinusoidal signal at the output of an 
oscillator never is truly monochromatic in the mathemat¬ 
ical sense, but has a line-width and noise side bands, has 
important implications. For example, it reduces the pos¬ 
sibility for a radio receiver to discriminate a week signal 
close to a strong one (Q)- 

The first paper on this subject I am aware of is chap¬ 
ter 15 in Edson’s “Vacuum Tube Oscillators” Q. In this 
work Edson assumed that an electronic oscillator could 
be considered an almost linear system whose nonlinear¬ 
ities gave a negligible contribution to the spectral com¬ 
position of the signal. 

Edson’s model for an oscillator was a sort of enhanced 
Q-multiplier brought to the limit of hugely amplifying 
and filtering the thermal noise always present in any real 
circuit. In this model the gain of the positive feedback 
loop was smaller than unity by an infinitesimal amount, 
to accomodate for the noise power at the ampliher in¬ 
put. The equivalent Q of the circuit resulted very very 
high, but not infinite, and its selectivity accounted for 
the Lorentzian shape of the generated signal spectrum, 
in agreement with experimental observations. 

This model could account for the presence of the nar¬ 
row noise skirt on sides of the generated signal (the car¬ 
rier) and for the 20 dB intensity decay per decade of 
frequency moving away from the carrier. 

This work did receive little attention and was little 
quoted in the literature in the years following. The au¬ 
thor himself declared that he “was never fully satisfied 
with the results from this approach ...” Q, to the point 
that, a few years after the book, he published a com¬ 


pletely new paper on the same argument [^. 

Moreover, experimental studies soon demonstrated the 
presence of a 1/f noise component in the frequency re¬ 
gion nearest to the carrier, that Edson’s model could not 
explain, not even including some amount of non linear¬ 
ity. Yet, Edson’s model had the great virtue of being 
very clear to understand and of directly connecting noise 
characteristics with electrical circuit parameters. 

Other works followed shortly after. In Leeson ap¬ 
proached the problem from an heuristic point of view, 
presenting an equation that described the commonly ob¬ 
served spectrum of feedback oscillators. This paper, 
though very popular in the following works, did not face 
the problem of connecting the observed noise spectrum 
properties with the electrical circuit characteristics. Lit¬ 
tle later Lax laid the basis for approaching the problem 
through statistics and interest into the argument was par¬ 
alleled by an analogous interest for the line width of a 
laser generated light beam Q. 

Recently interest into oscillator noise has greatly in¬ 
creased because of its relevance in the performance of 
telecommunication networks, both digital and analog, 
wired and wireless, as documented by the large number 
of recent contributions to the literature. See, e.g., 
and references therein. 

These works, as well as Lax’s one, greatly deepened the 
theoretical analysis of the stochastic processes behind the 
phase noise origin. They make large use of concepts and 
arguments from the world of statistics, at the expense, 
though, of somewhat missing a clearcut connection with 
the plwsical circuit, its parameters and its behaviour. 

In [9|, e.g., the analogy has been used between the 
oscillator signal phase fluctuations and the process of 
molecular diffusion. As suggestive and appropriate as 
it can be by principle such a connection, its application 
to the study of phase noise in oscillators can help only 
from a qualitative point of view. Otherwise, a double ef¬ 
fort would be required first to correlate circuit equations 
with molecular processes and then to bring results from 
molecular behaviour back into the circuits realm. 

A consequence of this state of affairs from a didac¬ 
tic point of view is that the subject is often relegated 
to books and courses of electronics at the highest spe- 
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cialization level or is proposed through an heuristic and 
qualitative only approach, following Leeson’s approach 
(see, e.g., [3). 

In this paper I will try to approach the problem making 
proht of the large and good amount of work already avail¬ 
able, through a simple but rigorous enough way, mak¬ 
ing use mostly of those concepts familiar to the world of 
electronics like Johnson’s noise, modulation, frequency 
conversion and vectorial (de)composition. Clearly, some 
mathematical extras to this menu are unavoidable, not to 
remain at the qualitative and descriptive level only. They 
are the Central Limit and the Wiener-Khinchin theorems 
and some properties of gaussian distributions. Luckily 
enough, currently they are well established and widely 
known theorems, also at the student level. 

The path I shall propose in this paper easily lends it¬ 
self to the implementation of numerical experiments, very 
useful from a didactic point of view to illustrate the the¬ 
oretical concepts and enhance confidence into the results 
obtained through equations. 

THE MODEL OSCILLATOR 

It is a common practice in a linear system to group all 
noise sources into a single one, conveniently placed and 
with an equivalent effect on the output signal. 

In this paper I have chosen not to follow this way for 
the sake of the simmetry both in the circuit and in the 
equations that will be written. The cost of this choice 
will be that most equations will contain the sum of two 
equivalent terms, one for current, one for voltage, with 
no further increase in complexity. The advantage will 
be a clearer splitting of noise into amplitude and phase 
components and the possibility to introduce into the cir¬ 
cuit some non linear behaviour without having to deal 
with the violation of the primary hypothesis “in a linear 
system” intrinsic to the common practice. In the last, as 
expected, it will be clear that the total noise power is the 
unique meaningful quantity in evaluating phase noise. 

In an electronic oscillator, like in every harmonic os¬ 
cillator, unavoidable sources of noise are the dissipative 
elements associated with the two energy storing compo¬ 
nents of the resonator, most often an inductance L and 
its series resistance R and a capacitance C and its parallel 
conductance G (fig. 1). 

The dissipative elements R and G are responsible both 
for the introduction of thermal noise into the circuit, a 
random voltage Vn{t) and current i„(t) associated with 
resistance R and conductance C?, and for the resonator 
damping. 

Actually, in a stationary oscillator values of damping 
and noise generation are disjoined from the physical val¬ 
ues of R and G, because of circuits and devices that 
provide the positive (regenerative) feedback indispens¬ 
able to add a negative damp that cancels out the natural 



FIG. 1: A typical resonator in an electronic oscillator is 
composed by an inductance L and a capacitance C, together 
with their resistance R and conductance G, that account for 
losses. The two generators in and are the noise sources 
inseparably linked to the dissipative elements. 

(positive) damp of the resonator and allows a persistent 
oscillation to be sustained. 

These circuits and devices introduce their own noise 
contributions that add to and/or Vn- 

As a consequence, the value of R and Vn and of G and 
in can be considered uncorrelated, only with the con¬ 
straint that the spectral density Vn{ui) and of Vn(t) 
and init) cannot be lower than what due respectively to 
the R and G of the physical circuit, namely 

Vn > \/AkTR In > \/VkTG 

where k is Boltzmann’s constant and T absolute temper¬ 
ature. 

So, it is possible to consider the practical case where 
damping is kept to zero by regeneration, i.e. R and G 
are nulled by equivalent negative terms, but some noise, 
both Vn and in, is present in the circuit. After putting 
some energy Wo into L and G there will be a persistent 
stationary oscillation. For the sake of simplicity, let’s 
consider the initial conditions v{t) = 0 1^ and i{t) = 1 A 
at t = 0 so that, in the absence of noise, oscillation would 
be 

v{t) = riosin(a;ot) 
i{t) = io cos{u!ot) 

with, as usual, wq = Stt/q = Ij'/LC and 

^Li'^it) + ^Gv'^it) = ^Lil = ^Gv^ = Wo 

In the state space described by the i,v circuit vari¬ 
ables in absence of noise the representative point S of 
the system describes, as well known, a circle with con¬ 
stant angular velocity Wq. In the presence of noise both 
the trajectory and angular velocity of the point are per¬ 
turbed, as described in fig. 2. 

Noise decomposition 

By effect of noise current in and voltage in a time 
dt charge upon capacitor G and current through inductor 
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FIG. 2: Decomposition of independent noise contributions 
dv = {in/C)dt and di = (vn/L)dt into the radial dP and 
tangential dT parts. Normalization of the i, v, di and dv 
variables over io and vo as expressed by eq. 13) is not shown to 
avoid cluttering the figure and should be implicitly assumed. 

L change by an amount (with effective G = 0 and R = 0) 
dv dt di = ^ dt (1) 

O J-j 

As Vn and in are statistically independent random vari¬ 
ables, so are dv and di. 

Whatever their instantaneous values, dv and di can 
always be decomposed into a tangential dT and a radial 
dP part as shown in fig. 2, with: 


holds. Thereafter, by transforming a couple of values 
(fn,*n) by eqs. 01, zeroing alternatively the dP or dT 
part and transforming back with the inverse of 01 , two 
couples of components (vntPnt) and {vnpPnp) are ob¬ 
tained that sum up to give the original values, 

but yield each one only a dT or dP perturbation: 


Vnt = —inRo sin P cos P + Vn COS^ p 

( 5 ) 

• 2 '^n . 

'^nt = sm if — — Sin (f COS ip 

Rq 

(6) 

Vnp = fn-Ro sin :/3COS:p-I-Un sin^ 

( 7 ) 

2 '^Tl • 

= Zn COS if — SlU (f COS if 

Ro 

(8) 

= x/A/G. 



Using these equations it is possible to decouple the dT 
and dP components of the random motion of the repre¬ 
sentative point S in the state space and go on studying 
their effect independently each other. 

Clearly, in a completely linear circuit as in fig. 1, such a 
decomposition is only a mathematical formalism. Actual 
voltage and current will be given only by the sum of the 
two components and will be the linear response of the 
LRCG network under the action of the in generators. 

But, when the circuit will be modified in order to sup¬ 
press one of the two noise components, namely the ampli¬ 
tude fluctuations, and this can be obtained only through 
the unavoidable introduction of a suitable non linear be¬ 
haviour, then the remaining tangential component (the 
phase fluctuations) will describe the true physical current 
and voltage. 


p2 



/ u \ ^ i'^L + v'^C 
21^0 


( 2 ) 


Phase fluctuations 


The effects of dT and dP over the motion of point S 
are independent each other: dP changes the oscillation 
amplitude P, dT its phase (p by the amount dp = dT/P. 

There is only a small link left between phase and ampli¬ 
tude fluctuations. If the noise generators in and Vn are 
stationary, as usual, their rms amplitudes remain con¬ 
stant. Hence, the same amount of noise from and 
yields a different effect on dp = dT/P, depending upon 
the instantaneous value of amplitude P. But this effect is 
vanishingly small indeed. In a real system the standard 
deviation of the P fluctuations is many orders of magni¬ 
tude smaller than the P average value: nobody would be 
concerned in a change of one part per million of a noise 
source relative intensity. 

By standard analytical geometry relations, the trans¬ 
formation 


By definition of angular frequency and using eq. 0 

, f in . Vn \ dt , , 

= —(9) 

This equation expresses how noise voltage and current 
from the and in sources induce fluctuations in the 
phase p of the oscillator signal. The two quantities vqCP 
and ifjLP represent the amplitude of the oscillation (the 
envelope) and its evolution over time. In any real os¬ 
cillator of practical interest, the amplitude is controlled 
by some kind of stabilization mechanism, as described 
in the next section; hence in eq. 0 amplitude drift can 
be neglected and P does not move from the initial value 
P = 1. This yields 

dp ,( in . Vn ^ 

w = — = Wo-k 1 • smyj-h — • COS(/3 I (10) 

The effect of the phase fluctuations is clearly the intro¬ 
duction of a frequency modulation over wq. In the con¬ 
ditions that are common in real (and useful) oscillators 
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this effect is very very small indeed, hence it is safe to 
write 


uJn = uj — ujQ = —^ • sin(u;oO + ■ cos{ujQt) (11) 

vqC iqL 

and 

(p{t) = UJot + ipn{t) ipnit) = f UJn{t')dt' (12) 

Jo 


(but in the numerical experiments the differential equa¬ 
tions of the circuit will be integrated without using such 
an approximation). 

If in and Vn are statistically independent sources of 
white noise, as it is the case for Johnson’s noise in re¬ 
sistors, they have a uniform spectral density. Hence, the 
product of their components with the sin() and cos() os¬ 
cillations at frequency wq has the effect of folding their 
constant spectral density over itself, yielding again a con¬ 
stant spectral density. Under such conditions, the spec¬ 
tral density as a function of frequency / of the resulting 
oscillation around frequency /o is the Lorentzian curve 
(normalized to unitary power over the full frequency 
range) 


^(/) 


1 

TTCr 


1 


f-fo 

a 


2 

+ 1 


(13) 


where a is the line-width of the generated signal. 

The path from eq. m to contains all the mathe¬ 
matical difficulties, so widely described in the literature. 
It is possible to travel this path in three steps, making 
use of two widely known theorems (Central Limit and 
Wiener-Khinchin theorem) and an established property 
of gaussian distributions. All of this will be used here 
without any demonstration, that is left to appropriate 
textbooks. 

The starting point is the requirement that in and Vn 
be white noise sources. It is not necessary they have a 
Gaussian distribution. Usually, uJn will not have a Gaus¬ 
sian distribution, apart some very special conditions, be¬ 
cause the product of in and with the sin() and cos() 
functions in eq. would alter their original Gaussian 
distribution, if any. In any case, ipn will have a Gaus¬ 
sian distribution, as stated by the Central Limit theo¬ 
rem, because ipn results from the sum of infinitely many 
independent random contributions, as seen in eq. (ini). 

The second step is to obtain the time autocorrelation 
function Rf{t) for the signal at the output of the oscil¬ 
lator, namely for the sin[v3(t)] or the cos[(^(t)] function. 

Following [nl [I 3 , using the complex notation, this can 
be written as 


Rf{t) = 


f 1 

—2- • Re < lim — 


^•Re 


T-j-oo T 






where Aq is the oscillation amplitude and the () brackets 
are a short for the {lim fj computation. 

Here the properties of the gaussian distribution play 
their part. Using the gaussian moment theorem, it is 
possible to write 




(14) 


A short but clear demonstration for this can be found 


12| and references therein. 


Evaluation of the exponent in the right hand side of eq. 
d is easier if it is made in a discrete time framework, 
like in the numerical integration of the circuit differential 
equations in the numerical experiments. It has been sup¬ 
posed that noise generators are stationary, hence time 
average can be replaced by variance (indicated by the 
overline): 


- <Pn(i + T)f) = [V5n(0) - ^n(T)f = (r) 

(15) 

At every At time step, ipn increases by the amount 
ujnJAt. Increments are random and uncorrelated, hence 
sum quadratically. After N = t/A d time steps 

= ujlAt T (16) 


It is to be recognized that 2uinAt is the spectral density 
of LOn (this will be made clear in the description of the 
numerical experiments): 


2ojlAt = 


1 

2 


v^oC^ 


1 

( _IL 

dWo VC 



(17) 


where has dimensions of s 

The final form of the autocorrelation function is 

= ^cos(a;or)e“^"^C (18) 

The third and last step is the application of Wiener- 
Khinchin theorem that states that the spectral distribu¬ 
tion F(u!) of a stationary random sequence is linked to 
the Fourier transform of its autocorrelation function: 


pCX) 

F{ui) = 4: / Rf{t) cos{u 
Jo 


r)dj 


Applying the theorem to function dT51) 


(19) 


pOC 

F{u;) = 2 Aq / cos(a;T) cos(wo''')f 

Jo 


-nlTi4 


= 


.-^lr/4 . 


{cos[(a;o — w)r] -I- cos[(aJo + uj)t]} dr 


(20) 
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The rapidly oscillating term with frequency wq gives a 
negligible contribution to the integral and can be dropped 
(the integral goes rapidly to zero when w moves away 
from wq)- Developing the integral of the remaining term, 
the Lorentzian of eq. (fT51) is obtained, with 


fj = 


D2 

Stt 


( 21 ) 


This result is illustrated in the Numerical Experiments 
paragraph, in the first numerical experiment, where only 
the components (u„t, int) from the full noise sources {vn, 
in) will be applied to the resonant circuit. The oscilla¬ 
tion amplitude will remain constant and the well known 
Lorentzian noise distribution will appear on sides of uiq. 


Amplitude fluctuations 

As anticipated above, every oscillator of practical inter¬ 
est has some kind of regulation mechanism that controls 
the oscillation amplitude. 

There is a twofold requirement for this. First, the con¬ 
dition for persistent oscillation is that the damping given 
by R and G be canceled out (it does not make a great dif¬ 
ference if each one separately, as considered in this work, 
or together as a whole). This is accomplished, as well 
known, using a positive feedback, that introduces into 
the circuit a negative R' and/or G' that exactly compen¬ 
sates the R and G damping effect. 

The word exactly is a term out of the realm of 
physics or engineering; however accurate a compensation 
is made, it will not be exact but for a short instant, there¬ 
after a drift of one kind or another will make damping 
move away from zero and make the oscillation rise or 
decay. 

Only a negative feedback regulation mechanism can 
continuously compare actual oscillation amplitude with 
a reference value and adjust the negative damping to con¬ 
trast any drift. 

Indeed, in a matematical sense damping can be made 
exactly null. This has been assumed in all the numeri¬ 
cal experiments, not inserting into the oscillator circuit 
either R or G. But in the second numerical experiment 
the two noise sources in and Vn will exchange energy 
with the LG resonator and make the oscillation ampli¬ 
tude fluctuate in a random way. 

The current generator for example, will induce a 
disturbance in the voltage across condenser G whose 
spectral amplitude will be 

14„(w) = Iniio) -(22) 

ojG -- 

(jJL 

At the resonance frequency ojq = Ij^/LG expression 
diverges, and indeed this is what would happen, 
given an infinitely long time. 


But the presence of the regulation mechanism, that 
cannot be avoided for the reasons explained above, keeps 
under control and limits the amplitude fluctuations also. 
And this is the second reason the regulation mechanism 
is needed. 

One of the simplest conceivable regulation mechanism 
is a first order control loop described by the equation 

^=K-{Po-P) (23) 

that has a stationary asynthotic solution P = Pq and a 
transient behaviour given by an exponential relaxation 
with time constant 1/K. 

The dynamical behaviour of the solution of this equa¬ 
tion is analogous to the oscillation amplitude decay in a 
resonant circuit with resonance frequency ujq and figure 
of merit Qq: 




with 



(24) 


Similar equations hold for io(t). 

Hence the behaviour we should expect for the noise 
induced amplitude fluctuations is the excitation of a res¬ 
onance around frequency ujg with a bandpass ±Aw = K. 
This process is akin to the virtual damping in Q and will 
be illustrated in the second numerical experiment. 

The scenario outlined up to now is completely consis¬ 
tent if it is recognized that considering only the (u„t, int) 
couple of components of the (z;„, in) noise in the circuit 
while discarding the (u„p, inp) components is fully equiv¬ 
alent to suppressing such fluctuations through a very ef¬ 
ficient amplitude regulation mechanism, as expressed by 
eq. (1^ . 


Upconversion of 1// noise 

From eqs. (unmii) it is clear that if in and Vn are 
coloured noise sources, also the resulting spectrum for 
LOn will be coloured, hence the main mathematical as¬ 
sumption in derivation of result (1131) will no more hold 
and the spectrum of the generated signal mostly will not 
be a Lorentzian curve. 

It is out of the purpose of this paper to analyze in any 
detail the implications of eqs. (IIOIIIII) with non-white 
and in noise sources. For a deeper analysis see, e.g., 0- 
Only an attempt has been made to illustrate the point, 
with the third numerical experiment where coloured noise 
sources have been used for Vn and in- 






6 


NUMERICAL EXPERIMENTS 

Numerical experiments have been carried out using the 
circuit of fig. 1 and numerically integrating the circuit dif¬ 
ferential equations in the presence of the two noise gener¬ 
ators in and u„. The computational engine of the Spice 
program is a very well optimized computer applica¬ 
tion to perform this task. In this work all computations 
have been made using version 24 of ngspice |14l |. an ac¬ 
tively maintained open source variety of Spice, and the 
jftw library [l^ to compute by Fast Fourier Transform 
(FFT) frequency spectra from temporal sequences (this 
library also available as open source). 

All computations that will be described can be afforded 
using a medium-high range PC, but their repetition to 
improve results quality by averaging noise spectra can be 
very time consuming. Hence, most of the repetitive work 
has been performed using an HP ProLiant DL560 Gen8 
machine equipped with 64 processors and 128 GB RAM. 


/at = l/(2At) = 50 kHz and the resulting noise spec¬ 
tral density is Vq/Jn = 2FjAt = 0.2 • 10“® V^ms/Hz 
(—67 dBV in a 1 Hz bandwidth). 

In a similar way the random variable of eg. 1111 with 
variance has spectral density = 2a;^At as given in 

eq.ini 

Generator bwn is a convenience to obtain an interme¬ 
diate result and can be used to draw distributions of Wn, 
as in fig. 5 below. 

Two further controlled generators (bvg and big) trans¬ 
form the independent voltages from vnv and vni into the 
pair of values Vnt, int given by equations (IM61) and play 
the role of and in in the circuit. 

Because the numerical values of L and C have been 
chosen to be equal, if the amplitude of the oscillation 
remains constant and equal to 1 (both V and A), with 
the given initial conditions sin ip and cos p can be written 
simply as 

sin p = i cos p = v 


The basic machinery used for the first numerical ex¬ 
periment can be seen in the following Spice list: 


*** experiment 1 *** 


11 0 2 
cl 1 0 

vnv 3 0 

vni 4 0 

bwn 5 0 

bvg 2 1 

big 1 0 


159.154943092e-6 ic=l 

159.154943092e-6 ic=0 

dc 0 trnoiseClOOm lOu 00000) 

dc 0 trnoiseClOOm lOu 00000) 

v=v(l)*v(3)-i(bvg)*v(4) 

v=v(l)*v(5) 

i=-i(bvg)*v(5) 


.control 
set nobreak 
set numdgt=12 

tran lOu 10000m Om 0.lu uic 
linearize 

print v(l) i(ll) > data.out 
quit 
. endc 


Values for L and C (11 and cl) have been chosen to 
give a resonance frequency /o as close as possible to 
1000 Hz io reduce the noise generated by truncation 
during FFT (in the experiments when phase is not af¬ 
fected) and have been set numerically equal to avoid un¬ 
needed computational overheads. This way equal ampli¬ 
tudes on vnv and vni give equal contribution to overall 
noise power. There should be no loss of generality in this 
choice. 

The two generators vnv and vni are white noise 
voltage sources that yield random sequences of values 
with a gaussian amplitude distribution and rms value 
Vg = 0.1 V. The cadence of the computed signals has 
been set to At = 10/is, hence Nyquist’s frequency is 


namely, sint/j = i(bvg) and cosp = v(l). 

First of all, the computational engine has been tested 
idle, i.e. performing computations with the amplitude of 
vnv and vni generators set to zero, for a time length of 
1 s. 



FIG. 3: Plot of i vs. v without noise (a, left) and with a 
noise of 0.1 Vrma {b, right). Computational length; 1 s (1000 
cycles). 


The results are shown in fig. 3 a, where values in the 
data file data, out have been plotted as i vs. v, i.e., i (11) 
vs. v(l). The plot covers 1000 cycles, each cycle con¬ 
taining 100 equally spaced samples and confirms phase 
coherence and amplitude stability of the computational 
engine over such a time length (indeed, a more accurate 
analysis shows a decay in amplitude of 0.04%/s). The 
noise floor due to the overall arithmetic accuracy is about 
— 145 dBV (or dBA) in a 1 Hz bandwidth. 


Experiment 1: phase fluctuations 

The computations described above have been repeated 
activating the noise generators vnv and vni at levels of 
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0.05, 0.1 and 0.2 Vrms over a time length up to 10 s, in or¬ 
der to obtain a frequency resolution of 0.1 Hz. The plot 
in fig. 3a changes, in that points do not repeat with the 
periodical cadence of 3.6°, rather scatter continuously all 
over the circumference. Neverthless the radius remains 
unitary, thanks to the behaviour of the bvg and big gen¬ 
erators that suppress the dP component from each couple 
of Vn, in values, leaving only the dT component (fig. 3&). 


> 

m 


E 

< 




100 mV 



freq [Hz] 


FIG. 4: Spectral amplitude in a 1 Hz bandwidth of the 
V signal across the C capacitor in the presence of the Vnt 
and int components only of the Vn and in noise generators. 
Noise generators intensity: 50, 100 and 200 mV. Results 
from the numerical experiment (dots) are plotted in compar¬ 
ison with Lorentz’s function (solid lines). The a parameter 
computed from eq. (I17II and (I21II is ct = 0.0785, 0.314 and 
1.257 Hz respectively from the lower to the upper curve). No 
time windowing has been applied with FFT, because the noise 
contribution due to truncation is 15-20 dB below the curves 
shown in the plot and becomes visible only for noise gener¬ 
ators amplitude of 25 mV or smaller. Each curve obtained 
averaging results from 15 repetitions of the experiment, with 
a time length of 10 s. 




FIG. 5: Distribution of the cUnAt values given by eq. (HID 
(on the left) and Bode plot of the right branch of the spectral 
amplitudes (on the right) obtained with two generators with 
uniform distribution (U) and with a single gaussian genera¬ 
tor (S). The third plot (G), for comparison, is made from the 
data of fig. 4. The expected analytical gaussian distribution 
is plotted as well together with the G curve, on the left, and 
is practically undistinguishable from the numerical distribu¬ 
tion. In all cases the total rms value of the generators was 
200 • V2 mV. 


The computations have been repeated again using for 
vnv and vni two generators with a uniform values dis¬ 
tribution instead of a Gaussian one, with rms amplitude 
200 mV (interval ±346.4 mV)\ 

vnv 3 0 dc 0 trrandom (1 lOu 0 346.4m) 
vni 4 0 dc 0 trrandom (1 lOu 0 346.4m) 

and once more with a single gaussian generator, either 
vnv or vni, with an rms amplitude 200 • mV. 

The distributions of the values given by eq. (fTTl) 
for the three cases are shown in fig. 5, on the left. Only 
the special case of the first computation exhibits a Gaus¬ 
sian distribution. Nevertheless, the other two cases also 
yield the same identical Lorentzian spectral amplitudes, 
as shown in the Bode plots on the right. In all cases, the 
Gaussian form of the (pn distribution is granted by effect 
of the Gentral Limit Theorem and can be seen examining 
the sum of a few consecutive values of UnAt. 


Experiment 2: amplitude fluctuations and 
stabilization 

In order to apply amplitude only fluctuations to the 
circuit, without affecting the phase, generators bvg and 
big are to be modified according to eq. (171181) . In this case 
P cannot be considered any more a constant, so that 


The spectral amplitude of the generated v signal across 
capacitor C is shown in fig. 4, together with the computed 
Lorentzian distributions, with the a values obtained from 
eq. (El) and (ED). The spectra for current i through 
inductance L are identical and are not shown. 


% V 

smp= — cos if = —- 

iqP vqP 

Thanks to the numerical equivalence of L and C the work 
can be done substituting the bwn, bvg and big generators 
in the Spice list of the first numerical experiment with 






bpq 

20 

0 

v=v(l)*v(l)+i(bvg)(bvg) 

bpp 

21 

0 

v=(v(4)*v(l)+v(3)*i(bvg))/v(20) 

bvg 

2 

1 

v=i(bvg)*v(21) 

big 

1 

0 

i=v(l)*v(21) 


Here also, generators bpq and bpp are only a conve¬ 
nience, to avoid the repetition of a relatively long com¬ 
putation common to both bvg and big. 

Fig. 6 shows the results obtained with such a ma¬ 
chinery with the P{t) plot (right) and the i vs. v plot 
(left). From the figure it is clear that now the phase 
proceeds at the constant angular speed of 3.6° per step 
(i.e., 1000.0 Hz) without any disturbances, while ampli¬ 
tude fluctuates by effect of energy exchange with the Vn 
and in generators. In the long time the oscillation am¬ 
plitude would fluctuate without a limit, as usual also for 
the phase, unless a phase lock mechanism to some time 
reference is used. 



FIG. 6: Amplitude fluctuations in the oscillation amplitude 
during 1 s with amplitude of vnv and vni noise generators set 
to 0.0114ms. 


In any case, an amplitude regulation mechanism is nec¬ 
essary for the reasons already explained. Even the purely 
computational approach of Spice applied to an ideal cir¬ 
cuit is not completely free from drift, as quoted in the 
introduction to the Numerical Experiments. 

The amplitude control expressed by eq. (1^ has been 
implemented computing the instantaneous deviation of 
energy in the circuit with respect to the reference value 
Wo 


and using this quantity as an error signal to actuate the 
radial only correction 

dP = KE{t)dt (25) 

This can be split into the dv and di components: 

dv = Vo cos • dP = v{t)dP 
di = iosimp-dP = i(t)dP 


Combining these expressions with equation o and (I25L 
the correction current and voltage 4 and Vc result: 

dv = ^dt = v(t)dP —>■ ic = KCv(t)£{t) (26) 

o 

di = ^dt = i(t)dP —>■ Vc= KLi(t)£{t) (27) 

Jj 

These corrections can be applied to the circuit adding 
their values to generators bvg and big, e.g. through gen¬ 
erator bpp, modifying its definition into 

bpp 21 0 v=(v(4)*v(l)+v(3)*i(bvg))/v(20)+ 

+ (sqrt(v(20))-l)*K 

and adding the command “. param K=. . . ” to set the 
appropriate value for the K constant. 

Computations have been performed with an amplitude 
of 10 mV for the vnv and vni generators and values of 
K from 10“^ to I. Results are shown in fig. 7, together 
with the expected resonance curves. The value of Qo 
for each value of K has been computed according to eq. 
(l24l) . taking into account that the C and L factors in eqs. 
(l26l) and (1271) actually are contained in the computational 
value of K: 


Qo 


LOo 

2K/C 


(28) 



FIG. 7: Spectral amplitudes in a 1 Hz bandwidth of the 
V signal in the presence of the Vnp and inp components only 
of the Vn and in noise generators. The noisy curves have 
been obtained with values of the K constant in the amplitude 
regulation mechanism of 1, 0.1, 0.01 and 0.001, from the lower 
to the npper curve (each one is the average of 25 repetitions 
of the numerical experiment, with a time length of 1 s). The 
smooth lines over the noisy ones are the resonance curves 
computed with the Qo values given by eq. (1281) . 


The frequency spectra clearly show the presence of a 
really monochromatic oscillation of amplitude 0 dBV, 
as set by the initial conditions, and the response of a 
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resonator with different values of Qo-, excited by the noise 
signals of and in- 

As already pointed out by most authors, an efficient 
amplitude regulation mechanism can reduce the contri¬ 
bution of the amplitude fluctuations to the overall noise 
to a level negligible in front of the phase fluctuations con¬ 
tribution. 

However, it is not correct to conclude that half the 
overall noise injected by and in into the circuit can 
be suppressed by the regulation mechanism. This hap¬ 
pens only in this very particular and abstract example in 
which the regulation mechanism uses two synchronized 
generators for voltage and current, that act only over the 
amplitude, avoiding any interference with the phase. 

While it is not impossible by principle to implement 
practically such a machinery it is not that easy. A pos¬ 
sible way, as suggested in (9|] is to make the regulating 
network, usually behaving as a voltage or a current gen¬ 
erator, to interact with the resonator only when the oscil¬ 
lation phase is at an appropriate point and the correction 
would require only a voltage or a current injection. 

Most often it happens that the regulating injected sig¬ 
nal contains a dT component over the dP one and this 
makes the noise contribution suppressed as an amplitude 
fluctuation to appear again as a phase fluctuation. 

Numerical experiments have been performed to verify 
this point (results not shown). Applying both the dT and 
dP noise components to the resonator with the amplitude 
control active (with K = 10), noise spectra identical to 
fig. 4 are obtained. Thereafter, suppressing the regulat¬ 
ing action of one of the bvg or big generators, correction 
is no more only radial and phase noise spectrum shows 
the resulting tangential increment. This is the common 
behaviour of real oscillators, where the splitting of the 
regulation signal into a dT and a dP components is not 
usually feasible. 


Experiment 3: upconversion of low frequency and 
1// noise 


Experiments in this section have been performed in the 
conditions of both experiment 1 and 2, i.e. with no am¬ 
plitude control, applying only the dT noise component, 
and with the amplitude control active (with K = 10), 
applying both the dT and dP components to the res¬ 
onator. This last condition has been implemented using 
the following definitions for the b. . generators: 


bpq 20 
bpp 21 
bvg 2 
big 1 


0 v=v(l)*v(l)+i(bvg)*i(bvg) 
0 v=(sqrt(v(20))-l)*K 

1 v=i(bvg)*v(21)+v(3) 

0 i=v(l)*v(21)+v(4) 


Identical results have been obtained with both imple¬ 
mentations of the experiment (apart a different level of 
the noise floor due to the arithmetics), because in both 


cases the dP perturbation component was effectively re¬ 
moved. 


Experiment 3a 

First, the oscillator behaviour has been studied under 
the application of a 10 mV low frequency sinusoidal sig¬ 
nal with frequency fi = 0.25, 0.5, 1, 2, 4, 8 Hz in place 
of the vnv and vni noise generators. With fi = 1 Hz the 
generators were 

vnv 3 0 dc 0 sin(0 10m 100) 
vni 3 0 dc 0 sin(0 10m 100) 

and the results are shown in fig. 8. 



freq [Hz] 


FIG. 8: Power spectral density of the v (and i) oscillator vari¬ 
able in the presence of a sinusoidal perturbation at frequency 
1 Hz and amplitude 10 mV for Vn and in- The computation 
has been extended over a time length of 20 s. A Blackmann- 
Harris window has been applied with FFT in order to reduce 
truncation effects. This spectrum has been obtained applying 
ony the dT perturbation component. The spectrum obtained 
the other way, i.e. with the amplitude control active, shows 
identical peaks and a noise floor 20 dB higher. 


Experiment 3b 


Thereafter, two coloured noise sources have been used 
for the Vn and generators, obtained summing the usual 
white gaussian noise with a 1//" noise component: 


vnvg 

3 

30 

dc 

0 

trnoise 

(0.1 lOu 0 

0 0 0 0) 

vnvf 

30 

0 

dc 

0 

trnoise 

(0 lOu 1.9 

Im 0 0 0) 

vnig 

4 

40 

dc 

0 

trnoise 

(0.1 lOu 0 

0 0 0 0) 

vnif 

40 

0 

dc 

0 

trnoise 

(0 lOu 1.9 

Im 0 0 0) 


A value close to 2 has been used for the n exponent in 
the 1//" distribution to improve the separation among 
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freq [Hz] 

FIG. 9: Power spectral density of the v (and i) oscillator 
variable in the presence of coloured noise sources in the role 
of the Vn and in generators. The plot inside the inset reports 
the spectral density of the noise generators (averages from 34 
computations 10 s long). 


the excess noise region and the oscillation frequency /o, 
while keeping an high enough excess noise amplitude. 

Results from this computations are shown in fig. 9, 
together with the Lorentzian curve associated with the 
white gaussian-only noise source. 

An analysis of results in fig. 8 and 9 will be given in 
the Discussion. 


DISCUSSION 


Results shown in fig. 4 and 7 are well known. Fig. 7 
shows that the fluctuations in the oscillation amplitude 
can be deeply reduced, to the level of being negligible, 
by a suitable circuit regulation mechanism. When the 
regulation is tight enough (AT > 1) only a flat noise floor 
is left, whose level decreases by 20 dB for each tenfold 
increase in K. 

The phase component of the noise induced signal fluc¬ 
tuations, instead, in the case of white noise sources, im¬ 
presses to the generated signal spectral distribution the 
shape of a Lorentzian curve (fig. 4). 

The most relevant and new aspect of this part of the 
work is the neat and direct mathematical path that starts 
from the splitting of the noise into the dP and dT compo¬ 
nents and leads to equations (1131), (EH) and (Oil) . These 
express both the functional form of the noise distribu¬ 
tion and the values of the constants it depends on, di¬ 
rectly expressed in terms of the unique two relevant cir¬ 
cuit characteristics: the energy Wq of the oscillation in 
the resonator and the rate of energy exchange (power) of 
the noise source(s) I^C + V^L. 

Remaining with the simple case considered in the first 


numerical experiment, in which 

In Vn 
vqC iqL 

eq. dm) becomes, e.g., 


n 


2 

n 


il/c 

2Wo 


(29) 


Substituting this last expression into eq. [HJ 

a = 

IGtt Wo 


(30) 


If the noise sources were only due to the losses in the 
resonator by effect of the parallel conductance G (or serial 
resistance i?, or both), one should have = AKTG and, 
by common formulas, Qo = uJoGjG. 

Combining this expressions into eq. dsni): 


27rv^GQo ^ ’ 

This result is identical to anal^ous results reported 
in the literature, e.g., eq. (25) in and eq. (70) in [^; 
but in the form of eq. (l30l) it is more pregnant and far- 
reaching. In a more general case, it is possible through 
eq. dm to sum up into the in and generators all noise 
contributions that exist in a circuit, then compute 
and the resulting line-width cr. From the details of the 
amplitude regulation mechanism, through eq. (131181) , it is 
possible to compute how much of the amplitude noise 
component is converted into phase noise, to be added to 
the native phase noise component. 

It is worth noting that the figure of merit Qo of the 
resonator does not appear in eq. dSOl), while it appears 
when passing to eq. EH). It has always been believed, 
and verified experimentally, that an high value for Qq is 
beneficial to reduce phase noise. Equations ([30)) and (|3Tll 
tell this is true, but only when the losses in the resonator 
are the main contribution as noise sources; otherwise, in 
passing from eq. (|30l) to (l3Tl) cannot be replaced by 
AKTG = 4:KTuJoG/Qo, but should become 


^2 ^2 , IKTujoG 

" Qo 


(32) 


where stands for all other sources of noise current dif¬ 
ferent from the parallel G of the resonator. This equation 
expresses the same concept as the Qioaded in 0. Also in 
Q it had been observed that an increase in Qo by any 
means, that also increase the noise, is vain. 

As already observed, what is truly relevant for the 
spectral purity is the noise power to signal energy ratio 
as expressed by eq. (EH- 


About the third experiment, the application of a small 
low frequency perturbation to the circuit, as made in 
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Experiment 3a, reveals the deeply non-linear behaviour 
induced by any efficient amplitude regulation mechanism. 
Either the radial component dP of the perturbation be 
removed a priori, as in Experiment 1, or be suppressed 
by the amplitude regulation, the identical result is the 
appearance of the same side-band components around 
the frequency of oscillation. 

The very interesting aspect to be observed is that this 
is not a simple effect of amplitude modulation (AM): the 
plot in fig. 9 clearly shows that lateral components are 
at frequency 1000 ± 2 Hz, ±4 Hz and ±6 Hz, hence 
at an offset frequency twice the 1 Hz frequency of the 
perturbation signal (and higher order harmonics). True 
AM sidebands can be obtained, instead, adding the mod¬ 
ulating signal to the expression sqrt(20)-l in the bpp 
generator definition (results not shown). 

A close analysis of the data in fig. 9 clearly shows that 
the central peak is no more exactly at 1000 Hz, as it 
sistematically was in all previous computations; rather, 
it is slightly moved to the left by a small amount that can 
be appreciated also in the plot. As a fact, the data.out 
file obtained from this experiment does contain a signal 
with average frequency of 999.95 Hz instead of 1000.00. 
This offset in the oscillator frequency depends upon the 
amplitude of the perturbing signal but does not depend 
upon its frequency. All frequencies from 0.25 to 8 Hz 
gave exactly the same offset; varying the amplitude in 
the 1 — 10 mV range a dependance of the offset with the 
square of the amplitude was observed. 

Reducing the value of the K constant while perform¬ 
ing the third experiment in the condition of Experiment 
2 , the sidebands and the frequency offset progressively 
reduce and vanish with AT —>■ 0. 

This makes sense, because with K ^ 0 every non¬ 
linear behaviour is removed from the circuit and only a 
linear resonator is left, excited by a voltage/current at 
frequency /;. In the signal spectrum only /o and fi can 
be found, with no mixing effect at all. 

The frequency doubling cannot be an artifact due to 
the particular amplitude regulation algorithm adopted; 
it appears identical with no regulation mechanism at all, 
only by effect of noise splitting and dP component sup¬ 
pression. 

These findings are to be kept in mind when observing 
the behaviour of an oscillator in the presence of a coloured 
noise source, as it is the case for 1// noise components. 
Results in fig. 9 (the noisy line) show the behaviour in 
such a condition, compared to the white noise behaviour 
of fig. 4. It is not a simple process of upconversion and 
appearance of an upper and lower band on side of /o; 
rather, it is a deformation and a shift to the left of the 
distribution and its peak due to a process of frequency 
pulling as described in Em. 

In view of these results, it is doubtful that this process 
can be considered the source (or the unique source) of 
1 // components observed in oscillators on side of the 


carrier. Other effects should be considered; for example, 
the direct modulation of the generated signal in passing 
through a component with equilibrium value affected by 
1// fluctuations [l7|. But this is argument for a new 
investigation and is out of the purposes of this paper. 
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